Can hydrogen be stored inside carbon nanotubes under pressure in gigapascal range? 
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By using a newly fitted multi-parameter potential to describe the van der Waals interaction 
between carbon and molecular hydrogen, we study the hydrogen storage inside carbon nanotubes 
(CNT's) under pressure in gigapascal range. Comparing with the results of graphite, we find that 
the shape change of the nanotubes (the curvature effect) provides a different storage mechanism for 
hydrogen. The negative free energy change for hydrogen storage inside CNT's makes it possible to 
use CNT's as the nanocontainer [Carbon 45, 315 (2007)]. 

PACS numbers: 61.46.-w, 61.48.-l-c, 62.50.-|-p, 64.70.Nd 



I. INTRODUCTION 

Technologies using hydrogen as an energy source are 
being developed rapidly, among which the hydrogen stor- 
age in carbon materials is one major aspcctji Dillon et 
al^ are the first that suggested the possibility of achiev- 
ing very high hydrogen storage capacity by using single- 
wall carbon nanotubes (SWCNT's). By experiments us- 
ing high-purity SWCNT's, hydrogen storage capacity of 
about 8 wt% at 80 K and 120 atm was reported by Ye 
et al-r- and about 4.2 wt% at room temperature and 100 
atm by Liu et al.A On the other hand, molecular dynam- 
ics simulations and the Monte Carlo method were used to 
study the hydrogen adsorptioni^i^ However, there are still 
many open questions to be resolved, for example, whether 
hydrogen is chemisorbed or physisorbed in carbon mate- 
rials and whether the intratube space of SWCNT's is 
more important than the intertube space. 

The report by Ma et alJ- that chemisorption is achieved 
by collision of H atom with energy of 1-3 eV, and the re- 
port by Chan et al^ on the breaking of the H-H bond 
between carbon nanotubes under high pressure show 
the possibility of the chemisorption. Very recently, an 
experimental study should be a breakthrough for the 
chemisorption)^ that carbon nanotube films are found to 
be hydrogenated and the C-H bonds can be broken by 
heating to 600° C. 

However, for the much more widely studied physisorp- 
tion, there are very few exciting results. In 2004, Chan et 
alJ^ found by ab initio calculation that H2 molecules can 
be well confined between graphene sheets under pressure 
in the gigapascal range due to the negative free energy 
change AG. This method is much better than many 
of the other simulations, because in those studies either 
the external pressure is not considered or the interac- 
tion between carbon atoms and H2 molecules are not 
well described. For example. Ma and Xia et aZi ^^'^^i^^ 
studied the energy between H2 molecules and the wall 
of the SWCNT and found helical structures of hydrogen 
inside the tube. However, we don't know whether the H2 
molecules want to stay inside the tube or not, because 
the tube in their studies is encapsulated. Furthermore 



their potentials describing the C-H2 and H2-H2 interac- 
tions are both Lennard-Jones (LJ) like, which are not 
accurate under high pressurei^ii^. 

In this study, by using molecular dynamics simulations 
with a more accurate multi-parameter potential for C- 
H2 van der Waals (vdW) interactioni^^ we investigate 
the hydrogen storage inside the graphene sheets and car- 
bon nanotubes under pressure. We find that the carbon 
nanotubes have a different storage mechanism from the 
graphite. For the graphite system, one must separate the 
graphene sheets by overcoming a large binding energy 
to store the H2 molecules. But in the nanotubes, such 
binding energy is easy to overcome because the storage 
of hydrogen relaxes the tube from fully collapsed shape 
(thin ellipse) to fat ellipse and provides more negative 
contribution to the total energy change. This is the main 
reason that hydrogen storage inside nanotubes might has 
a negative free energy change which is the key criterion 
for the storage. 

Furthermore, our research also shows the possibility of 
CNT-based nanocontainer in which the hydrogen can 
be filled inside the container under high pressure and can 
be locked inside the tube due to molecular valves located 
at the two ends of the tube. 



II. METHODOLOGY 

Energy changes - To calculate the free energy change 
we simulate pure molecular hydrogen (H2), pure car- 
bon system (C), and hydrogen-intercalated carbon sys- 
tem (C-I-H2), respectively. Thus the free energy change 
can be written as 

AG = AH - TAS = AE + pAV - TAS 

= AE + pAV - r(5(C+H2) - [5(C) + 5(H2)]) 
^ AE + pAV + TS{Ii2)- (1) 

Here we consider the entropy difference between pure and 
hydrogen-intercalated carbon system as zero, because un- 
der high pressure, the H2 molecules inside carbon systems 
arc well confined to form an ordered structure. 
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In order to study the hydrogen storage mechanism, we 
also divide the energy into four parts, 
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where E}^^ and E"^^ are the covalent and the vdW energy 
between carbon atoms respectively, -E'h2-H2 is the energy 
between H2 molecules, and -Eq-Hj is the energy between 
carbon and molecular hydrogen. Thus AE' can be written 
in the form 



l^E = A£;^c 



E, 



C-H, • 



(5) 



Systems - The cubic simulation cell for pure hydro- 
gen contains 2400 molecules. For hydrogen storage in- 
side graphite, we simulated one 1920-C graphite and 
7V-H2-intercalated 960-C graphite where N = 480, 960, 
and 1440, respectively. Periodic boundary conditions are 
used in all three directions (x, y, and z) for above simu- 
lation cells. 

For the tubes, one 600-C (15,15) SWCNT and one 400- 
C (10,10) SWCNT are chosen in our study. Inside the 
tubes, we put different number of hydrogen, such as 150, 
200, and 300 H2 inside the (15,15) tube and 80, 100, 
and 150 II2 inside the (10,10) tube. The nanotubes are 
periodic along their axial direction (z). 

Simulation details - The carbon atoms are coupled 
to a heat bath by using the Bcrcndscn velocity scaling 
method,— while all II2 molecules inside the carbon sys- 
tems are free. Simulations show the carbon atoms main- 
tain at the set temperature 300 K with reasonably small 
fluctuations and the II2 molecules can reach such tem- 
perature in picoseconds. For pure hydrogen, the temper- 
ature is controlled by the same method. 

For each direction with periodic boundary condition, 
the external pressure is applied using the Berendsen 
algorithmjii while for the xy plane of the carbon nan- 
otubes, a constant-pressure method for finite systems is 
usedji^ We also use another method to apply the hydro- 
static pressure on the tube in the xy plane, by introduc- 
ing quite a large number of inert gas (He) to transmit the 
pressure. These two methods produce the same results. 
However, the latter method is the real way to apply the 
pressure in experiments. 

The velocity verlct method is used to resolve the mo- 
tion equation, with the time step 1.0 fs. After the system 
reaches the equilibrium, all data are got by averaging in 
50 ps. 

Empirical potentials - At temperature of 300 K and 
moderate pressure range, H2 molecules can be treated as 
structureless spherical particles, and are modeled by the 
Silvera-Goldman potentialf^^ The C-C covalent bonds in 
graphite and carbon nanotubes are described by using 
the Brenner potential^ which has been widely used for 
carbon systems. For the C-C vdW interactions, a conven- 
tional treatment of LJ potential, V = 6*12/?'^^ — C^/r^, 



TABLE L Fitted parameters for the multi-parameter C-H2 
potential. 
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1166.92 eVA^" 
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where C12 = 2.48 x 10" eVA^^ and = 20 eVA^, is 
adopted. By using such potential, at 300 K and zero pres- 
sure, the binding energy for two graphene sheets is 30.7 
meV/atom with interplanar spacing of 3.30 A, which IS m 
agreement with both experimental and the density func- 
tional theory results those range from 25 to 57 meV/atom 
with the same spacing.— For the C-H2 vdW interactions, 
we use a more accurate multi-parameter potential^ 



(r) = exp(a — f3r — "fr^) 



C'lo , 

) ^^P 



r 



(6) 



which is fitted from ah initio results and can provide a 
good description for high pressure regimes. The param- 
eters for this potential are listed in Table |T1 



III. HYDROGEN INSIDE GRAPHITE 

Layered structure - We plot the hydrogen structure 
between two graphene sheets in Figure [TJ Layered struc- 
ture is formed under high pressure with the number of 
layer dependent on the molecule number. The top view 
[Figure [IJd)] shows all the H2 molecules are well con- 
fined between the hexagonal rings of the graphene sheets, 
in agreement with the Sandwich structure in Chan et 
al.^s studyiiS Figure [IJe) shows the trajectory of one H2 
molecule in the three-layered structure under 5.0 GPa. 
The molecule hops from one confined position to another 
in about 10 ps after the external pressure increases, and 
stays for ever. Different from the liquid state under the 
same pressure, the hydrogen is a solid because such hop- 
ping phenomenon is very few and can not be found any 
more under higher pressure. This is interesting because 
at room temperature the pure hydrogen becomes solid 
under a higher pressure about 8.9 GPai^ 

Simulations also show the hydrogen inside the graphite 
reaches this layered structure through two-step conden- 
sations. For example, in the graphite system with the 
number ratio of C:H2=2:1, the gas hydrogen undergoes 
the first condensation to form a two-layered structure 
under pressure about 0.5 GPa. With increasing the pres- 
sure, the hydrogen becomes more liquid-like. The second 
condensation happens under 2.0 GPa to form a one- 
layered solid-like structure. For the number ratio of 1:1 
and 2:3, the pressures for the second condensation are 
3.7 and 4.6 GPa, respectively. 
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FIG. 1: Hydrogen structure between graphene sheets under 6.0 GPa at 300 K. (a)-(c): different layers for different number of 
liydrogen; (d): a top view for the (a) structure; (e): the trajectory of one H2 molecule for the (c) structure under 5.0 GPa. 
The trajectory shows that the molecule hops from one confined position to another and stays there for ever. Such hopping 
phenomenon just happens under low pressure and at the time when the pressure increases. The length unit for the trajectory 
is angstrom (A). 



In Figure m we show the hydrogen volume (right) and 
the pressure of hydrogen along the xy plane (left). 

Energy changes - Table |ll] shows the energy changes 
calculated using Equation [T] and [5] The positive free en- 
ergy changes here are different from the ab initio study^i- 
mainly because of the different pAV values. Our studies 
show a negative pAV when the pressure is smaller than 
5.8 GPa, in agreement with the ab initio study. However, 
the values are different. For example, under 2.0 GPa the 
pAV term for C(H2)i/2 of -0.0146 eV per H2 molecule is 
much smaller than the ab initio result of -0.0975 eV. The 
reason is easy to find out that in the ab initio study, the 
supcrcell contains only 8 H2 molecules that even under 
small pressure the hydrogen is fully confined. 

We also provide the contributions to the total energy 
change AE to analysis the storage mechanism. One can 
find from the table that AE^q is nearly zero. This 



is because the graphene sheets keep their plane struc- 
ture all the time, without changing the covalent energy. 
The positive AE^q terms show that the storage process 
should overcome the binding energy between graphene 
sheets. Fortunately, the negative E'c-Ha is almost in 
the same magnitude as AEqq . The cancellation between 
these terms makes only Ai5H2-H2 left to the total energy 
change. A£'h2-H2 becomes more negative with increas- 
ing the pressure is because the hydrogen inside graphite 
is solid and can not be further compressed along the two 
directions {x and y) of the graphene plane due to the 
graphene's hexagonal structure, while the pure hydrogen 
can be compressed further in all three directions. 

In all three simulations, AE are negative beyond few 
GPa. However, the magnitudes of AE and pAV are so 
small that they can not overcome the quite large TAS 
term. In other words, such layered structure under high 
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TABLE II: Energy changes per hydrogen molecule for pressures under which hydrogen is fully confined between graphene 
sheets. All three simulations for different carbon-hydrogen number ratios show positive free energy change AG, which becomes 
more positive with increasing the number of hydrogen. Specially for the number ratio of C:H2=2:1, pAV, TAS, and the 
contributions to the total energy change AE are provided. 
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4.5 


0.0939 


0.1134 


0.1338 


0.0006 


0.1026 


-0.0279 


-0.0965 


-0.0211 


-0.0137 


0.1287 


5.0 


0.0820 


0.1050 


0.1118 


0.0010 


0.1000 


-0.0455 


-0.0903 


-0.0347 


-0.0088 


0.1256 


5.5 


0.0717 


0.0948 


0.1020 


0.0012 


0.0975 


-0.0617 


-0.0839 


-0.0468 


-0.0037 


0.1222 


6.0 


0.0619 


0.0858 


0.0935 


0.0014 


0.0949 


-0.0776 


-0.0769 


-0.0583 


0.0017 


0.1185 


6.5 


0.0522 
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0.0923 


-0.0936 


-0.0700 


-0.0695 


0.0071 


0.1146 


7.0 


0.0419 


0.0687 
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0.0021 


0.0895 
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-0.0629 


-0.0808 
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Pressure (GPa) Pressure (GPa) 

FIG. 2: The pressure (left) and the volume (right) of the 
hydrogen between graphene sheets. Green squares are taken 
from the experimental datai^i and the others are the results 
in our simulations. 

pressure is energetically favorable but not allowed due to 
the positive free energy change. 

Pressure of hydrogen - We also calculate the pressure 
of hydrogen between the graphene sheets, as shown in 
Figure [S] Though the pressure of hydrogen inside is 
smaller than the external pressure, more hydrogen in- 
tercalation is not allowed because of the positive AG. 
Negative AG is the first rule for hydrogen storage, then 
the pressure is another judgment that at the maximum 
storage the pressure of hydrogen is no larger than the 
external one. 



IV. HYDROGEN INSIDE CARBON 
NANOTUBES 

Hydrogen structure - Figure [3] shows the layered struc- 
ture of hydrogen inside the (15,15) CNT under high pres- 
sure. The hydrogen-intercalated system also undergoes 
step-by-step hydrogen condensations to reach the final 
layered hydrogen structure, similar to hydrogen inside 
the graphite. For the number ratio of C:H2=4:1, there is 
another phase transition under pressure about 6.1 GPa, 




6.0 GPa 6.2 GPa 

(al) C(H2)iy4 (a2) C(H2)iy4 (b) C(H2)iys (c) C(H2)iy2 (d) pure CNT 

FIG. 3: Hydrogen structure inside (15,15) CNT under 6.0 
GPa at 300 K. Specially for the number ratio of C:H2=4:1, 
there is a phase transition under 6.1 GPa, with different struc- 
tures before (6.0 GPa) and after (6.2 GPa) the transition 
shown in (al) and (a2). In order to show the curvature effect, 
the structure of pure CNT under 6.0 GPa is also shown. 



that before the transition only the middle part of the 
hydrogen is well confined while after the transition the 
hydrogen show one-layered structure in one part and two- 
layered in another. These structures are shown in Fig- 
ure [3l^al) and Figure [3l^a2), respectively. For the H2 
molecules, the main difference from the hydrogen inside 
graphite is the molecules at the two ends of the ellipse 
have more carbon neighbors and therefore have lower en- 
ergy than other H2 molecules in the middle of the ellipse. 

On the other hand and more important, with hydro- 
gen intercalation, the CNT wall relaxes itself to release 
its elastic energy. One can find the different bent shapes 
for the tube with and without hydrogen inside from Fig- 
ure [31 For the pure carbon CNT under pressure, the 
shape of the cross section is an elongated ellipse, with a 
large elastic energy at the two ends of the ellipse. But 
when the hydrogen is intercalated, the curvature becomes 
larger at the two ends. Such effect is new for CNT's, and 
will contribute more negative energy change for hydrogen 
storage. 

Energy changes - As we expected, hydrogen storage 
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Pressure (GPa) Pressure (GPa) 



FIG. 4: Energy changes for (15,15) CNT: C(H2)i/4. There 
are several jumps in the energy changes, corresponding to 
different phase transitions or hydrogen condensations: (i) the 
phase transition of pure (15,15) CNT (0.3 GPa), (ii) the first 
hydrogen condensation (0.9-1.1 GPa), (iii) the second hydro- 
gen condensation (2.2-2.4 GPa), and (iv) the final phase tran- 
sition of hydrogen to form a well confined structure (6.1 GPa), 
as described in the text. 



inside CNT's produces more negative AE. Fortunately, 
the magnitude of pAV is still quite small despite of the 
sign, thus the free energy change AG might be nega- 
tive under pressure. Figure |4] shows the energy changes 
for hydrogen inside (15,15) CNT with the number ratio 
of C:H2=4:li^ It is quite interesting to find that un- 
der pressure above 6.4 GPa, AG is negative (-0.0016 and 
-0.0101 eV per H2 molecule for 6.4 and 7.0 GPa respec- 
tively). For the number ratio of 3:1, AG is also negative 
under 7.0 GPa (-0.0011 eV). 

Considering that in the AG calculation the entropy 
difference between the pure and hydrogen intercalated 
carbon system is taken to zero, the real AG should be 
a little smaller, or in other words, more negative. This 
result is quite different from the graphite that under pres- 
sure in gigapascal range, hydrogen storage inside carbon 
nanotubc is free-energetically favorable. 

When more and more hydrogen (C:H2=2:1) is inside 
the tube, we find AG increases to 0.0187 eV per H2 
molecule even under 7.0 GPa. It shows the maximum 
storage of hydrogen in our pressure range is between the 
number ratio of C:H2=3:1 and 2:1. 

To find out the reason why hydrogen can be stored in 
CNT, we still should look into the contributions to the 
total energy change AE. In Figure [3] one can find that 
the covalent bonding energy between carbon atoms con- 
tributes now, shown in the AEqq term, as we have ex- 
plained above. Hydrogen inside the tube relaxes the tube 
wall and cause a negative energy change, which cancels 
the positive vdW binding energy change AE^^ . We will 
show that AE^f-, can be more negative when the tube 
radius decreases. However, for the left two contributions 
Ai?H2-H2 and i?c-H2J they are still as the same as those 
in graphite. In all, due to the curvature effect, hydro- 
gen storage inside CNT causes AE^^ to be negative and 



FIG. 5; Hydrogen pressure for the storage inside (15,15) CNT 
along the direction (z) of the tube axis (left) and the volume 
of the hydrogen, similar to Figure [2] 



makes the total energy change AE' large enough to over- 
come the TAS term. 

Pressure of hydrogen - In Figure [5] we show the cal- 
culation of the hydrogen pressure along the tube axis 
direction (z) and the volume for the hydrogen. We can 
also get the maximum storage capacity for (15,15) CNT 
from the figure. For the number ratio of C:H2=4:1 and 
3:1, AG is negative under ~ 7.0 GPa. At the same time, 
the pressure inside is smaller than the external one, thus 
more hydrogen can be stored due to the negative free 
energy change. 



TABLE III: AG for (10,10) and (15,15) CNT. x is the number 
ratio between molecular hydrogen and carbon atoms, labeled 
by C(H2)a;. The unit for AG is eV per H2 molecule. 

Pressure AG for (10,10) CNT AG for (15,15) CNT 



(GPa) a: = 1/5 1/4 3/8 x = 1/4 1/3 1/2 

0.6 0.0164 

0.7 -0.0063 0.0209 

0.8 -0.0204 0.0676 

0.9 -0.0414 -0.0048 

1.0 -0.0194 

3.8 0.0006 

3.9 0.0020 

4.0 -0.0001 

4.1 -0.0058 

6.3 0.0050 

6.4 -0.0016 

6.5 0.0051 

6.6 0.0004 

6.7 -0.0060 0.0005 0.0211 

6.8 -0.0067 0.0009 0.0204 

6.9 -0.0026 0.0022 0.0220 
7.0 -0.0101 -0.0011 0.0187 
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Smaller tube - When the radius of CNT decreases, 
the curvature plays more important role to the energy 
change. In our simulations for (10,10) CNT, we find that 
AG is negative even under <1.0 GPa, much smaller than 
the critical pressure for (15,15) CNT. Table Hill also shows 
that the critical pressure for a negative AG increases with 
increasing the number of hydrogen. 

V. CONCLUSION 

In summary, we have studied the hydrogen storage in- 
side graphite and carbon nanotubes. Three results are 
listed below: 

1) Hydrogen can not be stored inside graphite due to 
the positive energy change. The main reason is that 
the binding energy between graphene sheets is overcome 
only by the C-II2 vdW interaction and the H2-H2 en- 
ergy change. So the total energy change AE is not large 
enough to overcome the positive TAS term. The hydro- 
gen inside graphite is solid with hep structure, and has a 
smaller internal pressure than the external one. 

2) For the carbon nanotube, the hydrogen storage un- 
der pressure will relax the tube wall and bring an addi- 
tional contribution to the total energy change. So it is 
possible to achieve a negative free energy change. For 



example, under 7.0 GPa, hydrogen can be stored in- 
side (15,15) CNT due to the negative free energy change, 
with the storage capacity higher than 5.2 wt%. 

3) When the tube radius decreases, the curvature effect 
plays more important role in the total energy change. 
The negative free energy change can be achieved even 
under a small pressure. 

Our results also show the possibility to design the 
CNT-based nanocontainerii^ The hydrogen can be filled 
inside the tube under pressure in gigapascal range, and 
because of the molecular valves, when we unload the ex- 
ternal pressure, molecular hydrogen will be locked inside 
the tube. Such fiU-and-lock mechanism has been dis- 
cussed by Ye et oZ.,— and is proved applicable by current 
studies. 
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